Self-interaction effects in (Ga,Mn)As and (Ga,Mn)N 
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The electronic structures of Mn-doped zincblende GaAs and wurtzite GaN are calculated using 
both standard local-density functional theory (LSDA), and a novel pseudopotential self- interaction- 
corrected approach (pseudo-SIC), able to account for the effects of strong correlation. We find that, 
as expected, the self-interaction is not strong in (Ga,Mn)As, because the Fermi energy is crossed 
by weakly correlated As p - Mn d hybridized bands and the Mn 3d character is distributed through 
the whole valence band manifold. This result validates the extensive literature of LSDA studies 
on (Ga,Mn)As, including the conclusion that the ferromagnetism is hole-mediated. In contrast, 
the LSDA gives a qualitatively incorrect band structure for (Ga,Mn)N, which is characterized by 
localized Mn 3d bands with very strong self-interaction. Our pseudo-SIC calculations show a highly 
confined hole just above the Fermi energy in the majority band manifold. Such a band arrangement is 
consistent with (although by no means conclusive evidence for) a recent suggestion^^ that formation 
of Zhang- Rice magnetic polarons is responsible for hole- mediated ferromagnetism in (Ga,Mn)N. 

PACS numbers: Valid Pacs 



I. INTRODUCTION 

The discovery of ferromagnetism in Mn-doped GaAs 
and InAgi has changed conventional thinking about elec- 
tronics. In fact these diluted magnetic semiconductors 
(DMSs) combine the functionalities of semiconductors 
with those of magnetic materials, paving the way for a 
complete integration of data storage and logic on the 
same electronic componenli^. Unfortunately the most 
widely studied DMS, (Ga,Mn)As, shows a Curie temper- 
ature {Tc) far below room temperature and therefore ap- 
pears to be unsuitable for commercial applications. This 
has led to a considerable world-wide effort in searching 
for new DMSs with higher T^s. 

Part of this research is guided by a pioneering theo- 
retical work by Dietl and collaborators, who predicted 
remarkably high T^s for wide-gap magnetic semiconduc- 
tors at moderate Mn doping^. Their model is based on 
the Zener mechanism for ferromagnetism, according to 
which a conduction electron can hop across two magnetic 
atoms through a double exchange of place with the va- 
lence electrons of the bridging anions. Spin-conservation 
implies that this conduction can only occur between spin- 
aligned magnetic centers, thus a ferromagnetic ordering 
of the Mn ions is expected. In the model, Tc is crucially 
dependent on hole and Mn concentration, and on the 
strength of their mutual coupling. High Tc ferromag- 
netism in wide-gap DMSs has recently been confirmed 
by several groups who reported the existence of above- 
room-temperature ferromagnetism in Gai_2;Mna;N^i^iSii. 
However, in contrast to "conventional" As-based DMSs, 
where the existence of double-exchange induced ferro- 



magnetism is clear, the origin and the nature of the 
observed magnetism in Gai-j^Mn^jN are controversial. 
On one hand, both p- and n-doped samples are, appar- 
ently, ferromagnetic, thus indicating that hole-mediated 
double-exchange may not be, in fact, the driving mecha- 
nism. More importantly, recent magnetic circular dichro- 
ism (MCD) measurements^ on high Tc samples show 
that the MCD signal in the vicinity of the GaN band- 
gap comes from a paramagnetic phase. This suggests the 
possibility of two phases in (Ga,Mn)N: a paramagnetic 
phase, attributed to substitutional Mn at the Ga site, 
and dominated by s,p-d interactions between Mn and 
the host GaN, and an unidentified ferromagnetic phase, 
not detected by MCD, which might be due to small par- 
ticles that are invisible to diffraction methodsS. 

In summary, there are sound indications that 
(Ga,Mn)As and (Ga,Mn)N may not be, in fact, very 
similar as far as magnetic ordering is concerned. Thus 
a careful examination of their microscopic properties is 
necessary in order to give a correct interpretation of the 
experimental results. 

With this goal in mind, in this paper we investigate 
the electronic structures of zincblende (Ga,Mn)As and 
wurtzite (Ga,Mn)N using first-principles density func- 
tional theory (DFT) calculations. Band energy calcu- 
lations for these compounds obtained within the familiar 
local-spin density functional theory (LSDA) are already 
present in the literatureiiSiiiiiSii^ Here we present the 
results obtained with our recently developed pseudopo- 
tential self-interaction free approach (pseudo-SIC).^ 

It is known that the band energies calculated within 
LSDA have systematic errors caused by the presence in 
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the one-electron Kohn-Sham (KS) potential of its self- 
interaction (SI)^^ (that is the interaction of an electron 
with the potential generated by the electron itself). The 
SI is especially relevant in ionic and mixed ionic-covalent 
materials. Indeed the poor LSDA description of the elec- 
tronic properties of II- VI semiconductors^^, transition 
metal oxidesi^, and magnetic perovskitesiSi is a conse- 
quence of the SI^i. In general, any compound character- 
ized by localized electron charges is, in principle, affected 
by a non- vanishing, unphysical, SI contribution. 



II. METHODOLOGY 

In this paper we use our recently developed self- 
interaction free local-density scheme^ based on ultrasoft 
pseudopotentialsi^ and a plane wave basis. The scheme 
is described in detail in Ref-A; here we summarize the key 
points. The Kohn-Sham equations are modified in order 
to effectively extract from the one-electron potential the 
spurious self-interaction contribution: 



Regarding the compounds investigated in this paper, 
it is known that the LSDA band structure of bulk GaN 
shows severe discrepancies with photoemission results 
(for example the fundamental energy gap is severely un- 
derestimated and the position of the d bands is incor- 
rect). Furthermore, the LSDA is generally inaccurate 
in describing the electronic properties of oxides of man- 
ganese, and in particular fails to give the correct position 
of the Mn d states with respect to the oxygen-derived 
valence manifold. 

Since the band gap of the host material and the posi- 
tion of the transition metal d states are crucial in the de- 
scription of these diluted magnetic semiconductors, there 
is the necessity to confront the LSDA description with al- 
ternative approaches. Our pseudo-SIC is ideally suited 
for this task since it repairs to a large extent the LSDA 
failures and provides a reliable description of the elec- 
tronic properties for a large range of materials, including 
GaN and Mn-based perovskites. In addition, and in con- 
trast with other beyond-LSDA approaches, our method 
preserves the conceptual simplicity and the computa- 
tional feasibility of the LSDA, and thus the possibility 
to approach large-sized systems such as the diluted mag- 
netic semiconductors. 

In this paper we determine the extent to which correc- 
tion of the SI changes the electronic properties and the 
chemical pictures of wurtzite (Ga,Mn)N and zincblendc 
(Ga,Mn)As. We find that in (Ga,Mn)As the localization 
of the Mn 3d electrons is weak because both on-site and 
off-site hybridizations occur, and, as a consequence, the 
SI is not strong. Therefore the LSDA already provides 
a good description of (Ga,Mn)As and the pseudo-SIC 
does not make large qualitative changes. In (Ga,Mn)N, 
however, the SI is much stronger, due to the large N 
electronegativity, and consequent higher ionicity of GaN. 
Therefore the LSDA description is inadequate and the 
pseudo-SIC gives qualitatively different results. The dif- 
ference in localization between the two compounds in 
turn has a profound effect on their electronic properties. 

This paper is organized as follows: in section ^ we 
describe the methodology and the technicalities used for 
the calculations. Sections IIIII and Hvl present the results 
for (Ga,Mn)As, and (Ga,Mn)N, respectively. Finally, in 
sections IVI and IVIl we present our conclusions. 
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where Vpp is the usual ion-core pseudopotential projec- 
tor, V-^^^ the LSDA-KS potential, and V^^^ the SIC term. 
The latter is cast as a nonlocal Kleinman-Bylander-type 
pseudopotential projector: 
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Here (j)i{r) is the pseudowavefunction, ^•^^^.i^'i i''^)] the 
SI potential of the j*'' atomic orbital {i runs over angu- 
lar quantum numbers and atomic positions, cr over the 
spin), and p"^ is the corresponding occupation number 
calculated self-consistently by projecting the Bloch state 
manifold onto the basis of atomic pseudowavefunctions. 
This procedure is no more computationally demanding 
than the usual LSDA, and, in contrast with other SIC 
schemesi^, it can be applied to both metals and insula- 
tors, a feature which is crucial for the present application. 

In earlier work our method was applied to a number of 
different systems including III-V and II- VI ionic semicon- 
ductors (e.g. GaN and ZnO),-^ transition metal monox- 
ides (MnO, NiO)^ and perovskites (CaMnOs, BaTiOa, 
YMnOa)^^. In all cases the pseudo-SIC obtained both 
the correct physical structures, and band structures in ex- 
cellent agreement with photoemission experiments. This 
series of successes lends strong support to the credibility 
of the present results. 

We use ultrasoft pseudopotentials to allow moderate 
cut-off energies (35 Ryd and 40 Ryd for (Ga,Mn)As 
and (Ga,Mn)N, respectively). For the pseudo-SIC cal- 
culations the Ga 3d electrons are treated as core states 
(this choice is fully justified in Ref.^.) For zincblende 
(Ga,Mn)As and wurtzite (Ga,Mn)N 16-atom and 32- 
atom supercells are used respectively, in each case with 
one Mn ion substituting for one Ga ion, corresponding to 
doping levels of ~ 12% and ~ 6%. The lattice parameters 
are fixed to the values calculated for the host materials 
(a = 10.68 bohr for GaAs, a = 6.03 bohr, c= 9.802 bohr, 
u= 0.377 for GaN), which are in good agreement with 
the experimental data. 



3 



> 



u 

e 



GaMnAs LSDA 
majority minority 
3 
2 
1 

-1 
-2 
-3 
-4 
-5 
-6 



M 



-7 
-8 



M 



FIG. 1: Band structure of zincblende (Ga,Mn)As calculated 
within LSDA for majority (left) and minority (right) spins. 



III. ZINCBLENDE (GA,MN)AS 
A. LSDA 

The calculated LSDA band structure of (Ga,Mn)As 
(in figure^ is that of an half-metal with metallicity in 
the majority spin band. Here the most important fea- 
ture for the electronic and magnetic behavior is given by 
the three bands crossing the Fermi energy in the major- 
ity manifold. They are rather dispersed and span ~2 
eV between F and M. An orbital decomposition of these 
bands shows that they consist of Mn d ^2 character, as 
well as Mn and As p, with significant contributions from 
both nearest and next-nearest neighbor As ions. This 
indicates that part of the charge is extended beyond the 
Mn-centered tetrahedron. The next set of bands lower in 
energy corresponds to the valence band top of bulk GaAs 
at the F point {E — — 0.73eV) but mixes in Mn charac- 
ter elsewhere in the Brillouin Zone, as does the doublet 
below it. An orbital decomposition shows that the Mn d 
character is not restricted to a single set of bands close to 
the Fermi level, but instead is widely distributed through 
the GaAs sp valence manifold. 

In the upper part of Table we report the occupa- 
tion numbers of individual orbitals. The total magneti- 
zation per cell is '~ 4 /^s, almost entirely due to the Mn 
(3.72 hb)- The Mn d electron charge is 5.48, of which 
4.61 and 0.87 comes from spin-up and spin-down elec- 
trons, respectively. 

In conclusion, according to the LSDA description (and 
in agreement with previous calculations^"'") (Ga,Mn)As 
is dominated by Mn d - As p hybridization rather than 
on-site electron correlation. However, since the overesti- 
mation of p-d hybridization and the suppression of on-site 
Coulomb energy are general tendencies of the LSDA, a 
comparison with the pseudo-SIC is necessary in order to 



TABLE I: Occupation numbers for selected orbitals of 
(Ga,Mn)As calculated within LSDA (top) and pseudo-SIC 
(bottom). The arrows indicate spin- up (majority) and spin- 
down (minority) components, respectively. Asj^jj indicates 
the As atoms adjacent to the Mn ion. 
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FIG. 2: Band structure of zincblende (Ga,Mn)As calculated 
within pseudo-SIC for majority (left) and minority (right) 
spins. 



assess the reliability of this picture. 



B. pseudo-SIC 

In Figure El we show the band structure calculated 
within pseudo-SIC. It is immediately apparent that the 
main features are not changed significantly with respect 
to the LSDA, apart for the expected increase of the fun- 
damental energy gap of the bulk GaAs due to the down- 
ward shift of the occupied valence band manifold with 
respect to the empty bands. Also the band manifold 
occupies an energy region ~ 1 eV wider than that in 
the LSDA result. Interestingly, the three majority, half- 
occupied bands have an even larger dispersion than in 
LSDA, but the shape of the bands is very similar. 

In the lower part of TableQlwe report the orbital occu- 



4 



pations calculated within pseudo-SIC. The total magne- 
tization is 4 iiB, as in LSDA. However, as a consequence 
of a larger Mn d spin splitting, the Mn magnetic mo- 
ment (4.58 iib) is larger than in LSDA, but it is com- 
pensated by a non discardable spin-polarization on the 
As (0.17 ^b) which is antiparallel to the Mn moment. 
(Table^reports the values for a single orbital (p^,) which 
has a magnetization of -0.06 fiB in pseudo-SIC. The total 
magnetization for p^: -I- Py -I- Pz is -0.17 /is.) 

In conclusion, the pseudo-SIC result is a sound con- 
firmation that the p-d hybridization and the weakly- 
correlated character of the Mn d electrons is in fact a 
reliable prediction of the LSDA, and validates the exten- 
sive earlier computational literature on (Ga,Mn)As. (for 
a review see RefiiS). 



> 

u 
e 



majority 



LSDA 



minority 




IV. WURTZITE (GA,MN)N 
A. LSDA 

In the hexagonal crystal symmetry of the wurtzite 
structure, the d orbitals are split into two doublets {{dxy, 
dx2-y2), and {d^z, dyz)) and one singlet (d^), where x 
and y lie in the hexagonal plane and z is parallel to the 
c axis. Similarly, the p orbitals are split into a doublet 
{px, Py) and a singlet (pz). In Figure 13 we report the 
band structure of wurtzite (Ga,Mn)N calculated within 
the LSDA. Like (Ga,Mn)As, this material is half-metallic 
and has three partially filled bands close to the Fermi 
energy. However this three-band manifold (consisting 
of a doublet composed of Mn dxy-dx2-y2 and the p^:- 
Pj, orbitals of the three planar N, and a singlet derived 
from Mn d^^-pz and on-top N pz) is much narrower in 
(Ga,Mn)N than in (Ga,Mn)As. The doublet below these 
states {E = -1.39 eV) at F is also a p-d combination 
(of Mn dxy-dx2-y2^ and N p^;) with a strongly localized 
electron charge confined almost entirely to orbitals in- 
ternal to the Mn-centered tetrahedron. The small band 
dispersion (less than 1 eV) is evidence of the strongly- 
correlated character. The same features can be appre- 
ciated in the orbital-resolved density of states, reported 
in Figure^ (only majority spin components are shown). 
The hybridization between Mn d and the p orbitals of 
the four surrounding N is clearly visible, while almost no 
contribution comes from second neighbor N. 

The orbital occupation numbers for (Ga,Mn)N are 
given in the upper part of Table The Mn magnetic 
moment is 3.63 fiB with a total Mn d electron charge 
of 5.3. Notice how the much stronger electronegativity 
of N with respect to As is reflected in the comparison 
between the orbital occupations in (Ga,Mn)As (Table P) 
where roughly 60% and 30% of the s, p electron charge is 
distributed to As and Ga, and in (Ga,Mn)N, where the 
proportion is 70% to N and 20% to Ga. 

These band structures agree qualitatively with earlier 
LSDA calculationsSSiSi, although the exact position of 
the Mn d levels differs according to the details of the 



FIG. 3: Band structure of wurtzite (Ga,Mn)N calculated 
within LSDA for majority (left) and minority (right) spins. 
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FIG. 4: Orbital-resolved density of states of wurtzite 
(Ga,Mn)N calculated within LSDA. Only majority-spin or- 
bitals are shown. For Mn, light-shaded and dark-shaded ar- 
eas represent the contributions from doublet (da;z,dt,z) and 
singlet d^2, respectively. 



pseudopotentials or the lattice parameters used in the 
calculations. 



B. pseudo-SIC 

Strong electron localization makes the LSDA re- 
sults for (Ga,Mn)N potentially affected by a large self- 
interaction. Our calculations indeed confirm these expec- 
tations. The pseudo-SIC band structure and the orbital- 
resolved density of states are shown in Figures [S] and |^ , 
respectively. 

Due to the self-interaction, the fully-occupied dxz, and 
dyz bands are pushed down in energy, well below the 
p-like valence band top of the host material. Around 
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FIG. 5: Band structure of wurtzite (Ga,Mn)N calculated 
within pseudo-SIC for majority (left) and minority (right) 
spins. 
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FIG. 6: Orbital-resolved density of states of wurtzite 
(Ga,Mn)N calculated within pseudo-SIC. Only majority-spin 
orbitals are shown. 



the Fermi energy two electrons occupy the other three 
majority Mn d bands, similarly to what we have seen 
within LSDA. However now these three bands are closer 
to the top of the p band manifold, and, due to the large 
increase of the GaN energy gap, much farther from the 
bottom of the s-like conduction bands. 

In the lower part of Table ^ the occupation numbers 
are shown. The local magnetic moments are 4.18 /is for 
Mn, 0.06 hb for the on-top N, and -0.10 for each of 
the three planar N. This means that the coupling between 
the on-top N pz orbital and the majority d^2 orbital is 
ferromagnetic, whereas the coupling between planar N px 
and Py orbitals with dxy and da;2_j,2 orbitals is antipar- 
ralel, thus contributing negatively to the total magnetic 
moment M=4 ^b- 

This remarkable increase of Mn magnetic moment with 



respect to the LSDA values is a typical self-interaction 
correction effect: the spatial localization is increased and 
the hybridization reduced. 

Finally, notice that these results are not very sensi- 
tive to the Mn concentration: band energy calculations 
of wurtzite (Ga,Mn)N with 12% Mn doping within both 
LSDA and pseudo-SIC shows no relevant changes with 
respect to the corresponding calculations at 6% concen- 
tration. 



V. DISCUSSION 

The chemical picture of (Ga,Mn)N obtained within 
the pseudo-SIC (which we beheve to be the most accu- 
rate representation) is markedly different both from the 
LSDA description of (Ga,Mn)N, and also the behavior 
of (Ga,Mn)As. The difference between (Ga,Mn)N and 
(Ga,Mn)As can be traced back to the large electroneg- 
ativity of N, which in turn favors ionicity and charge 
transfer over hybridization. This aspect, already evident 
within LSDA, is further emphasized within pseudo-SIC, 
due to its tendency to favor electronic configurations with 
completely filled or empty orbitals>22i 

The differences are so large that they will certainly 
lead to different charge mobilities in the two materi- 
als, and affect the nature of the magnetic ordering as 
well. Indeed the band alignment obtained within pseudo- 
SIC for (Ga,Mn)N may not be appropriate for validat- 
ing the usual Zener picture for ferromagnetism. In fact, 
according to our calculation, there are no free holes in 
the GaN valence band. Instead, the holes are local- 
ized near the Mn impurity, and occupy a triplet of 1 eV 
wide bands around the Fermi energy. This could lead to 
two possible scenarios. On the one hand substitutional 
Mn in GaN could be responsible only for the paramag- 
netic component of the magnetization as indicated by 
MCD measurements^. On the other hand this could be 
indicative of the formation of a Zhang-Rice polaron in 
(Ga,Mn)N as has been recently suggested-^. According 
to this model, Mn doping in GaN behaves as an effective 
mass acceptor, with a configuration of d^ plus a hole lo- 
calized on the magnetic ion. This local singlet could move 
through the sublattice of Mn^+ ions and cause ferromag- 
netic ordering (still mediated by a double-exchange mech- 
anism) at much higher Tc then that of the usual Zener 
"free hole" double exchange. Our calculations show a 
magnetic moment for Mn of 4.18 /x^, which is sugges- 
tive of a d^ configuration. However the presence of very 
flat, partially occupied bands with d character allow us 
to re-interpret our results in terms of a Mn d^ configura- 
tion plus a tightly bound hole of d-character. These are 
the conditions for the formation of a Zhang-Rice polaron. 
To further assess the validity of these hypotheses, larger- 
sized calculations, capable of describing the properties 
of these compounds in different magnetic orderings (e.g. 
paramagnetic and antiferromagnetic) will be necessary in 
the future. 
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TABLE II: Occupation numbers for selected orbitals of 
wurtzite (Ga,Mn)N calculated within LSDA (top) and 
pseudo-SIC (bottom). Only the 4 nitrogens surrounding Mn 
are considered: N^Qp is that on-top of Mn, N^-^^gg s one of 
the three N at the vertices of the triangular base. The arrows 
indicate spin-up (majority) and spin-down (minority) contri- 
butions, respectively. 

Ntop Nbase Mn 

T i T i T i T I 

p^ 0.24 0.24 0.73 0.71 0.70 0.70 0.20 0.15 

p^ 0.23 0.23 0.72 0.71 0.74 0.71 0.20 0.15 

d^y 0.85 0.20 

d^^ 0.93 0.14 

d^2 0.82 0.24 

p^ 0.22 0.22 0.77 0.75 0.67 0.72 0.20 0.17 

p^ 0.22 0.22 0.79 0.75 0.74 0.74 0.20 0.15 

d^y 0.88 0.12 

d^^ 0.95 0.08 

d^2 0.98 0.16 



VI. CONCLUSIONS 

In summary, we have compared the electronic structure 
of ferromagnetic (Ga,Mn)As and (Ga,Mn)N within stan- 



dard LSDA and within our pseudo-SIC approach which 
sets the KS potential free of the spurious self-interaction 
part. Our results indicate two main conclusions. First, 
that the LSDA picture of (Ga,Mn)As as a half-metal, 
with the Fermi energy crossing three weakly correlated 
bands which derive from Mn d - As p hybridized states, 
is indeed a good description, largely unchanged when the 
self-interaction is corrected. However for (Ga,Mn)N the 
LSDA description is inadequate because the d electrons 
are strongly correlated, and thus the pseudo-SIC descrip- 
tion is markedly different from that of the LSDA. Second, 
that the behavior of (Ga,Mn)N is qualitatively differ- 
ent from that of (Ga,Mn)As. Three flat bands of d^^j,, 
da;2_j,2, and d^^ character are partially occupied by two 
spin-polarized electrons, confined within the Mn-centered 
tetrahedron, thus effectively describing a highly confined 
hole around the Fermi energy which may be more consis- 
tent with the formation of a Zhang-Rice polaron, rather 
than with the free-hole mediated double-exchange mech- 
anism of the Zener model. 
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